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ABSTRACT 

We study the X-ray properties of 393 optically selected early-type galaxies (ETGs) over 
the redshift range of z ~ 0.0-1 ,2 in the Chandra Deep Fields. To measure the average 
X-ray properties of the ETG population, we use X-ray stacking analyses with a subset 
of 158 passive ETGs (148 of which were individually undetected in X-ray). This ETG 
subset was constructed to span the redshift ranges of z = 0.1-1.2 in the «4 Ms CDF-S 
and «2 Ms CDF-N and z = 0.1-0.6 in the w250 ks E-CDF-S where the contribution 
from individually undetected AGNs is expected to be negligible in our stacking. We find 
that 55 of the ETGs are detected individually in the X-rays, and 12 of these galaxies 
have properties consistent with being passive hot-gas dominated systems (i.e., systems 
not dominated by an X-ray bright Active Galactic Nucleus; AGN). On the basis of 
our analyses, we find little evolution in the mean 0.5-2 keV to 5-band luminosity ratio 
{Lx/ L b [1 -fz]*'2) since z « 1.2, implying that some heating mechanism prevents the gas 

from cooling in these systems. We consider that feedback from radio-mode AGN activity 
could be responsible for heating the gas. We select radio AGNs in the ETG population 
using their far-infrared/radio flux ratio. Our radio observations allow us to constrain the 
duty cycle history of radio AGN activity in our ETG sample. We estimate that if scaling 
relations between radio and mechanical power hold out to z w 1.2 for the ETG population 
being studied here, the average mechanical power from AGN activity is a factor of «1.4— 2.6 
times larger than the average radiative cooling power from hot gas over the redshift range 
z « 0-1 .2. The excess of inferred AGN mechanical power from these ETGs is consistent with 
that found in the local Universe for similar types of galaxies. 
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1 INTRODUCTION 

The most successful theoretical models of galaxy evolution (e.g. 
Bower et al. 2006 and Croton et al. 2006) require that feedback, in 
the form of energetic outflows from active galactic nuclei (AGNs), 
will have a fundamental influence on the evolution of intermedi- 
ate and massive galaxies. In these models, energy injected from 
AGN radio Jets heats the interstellar and intergalactic mediums 
of massive early-type galaxies (ETGs) and further drives inter- 


stellar gas out of these systems. This energy injection effectively 
quenches star formation and supermassive black hole (SMBH) ac- 
cretion and prevents galaxies and SMBHs growing. The interstel- 
lar gas itself is thought to be produced by evolving stars eject- 
ing material through stellar winds and supemovae (at a rate of 
« 1.3[Z,5/10*’i5_3]Afe>r‘'*; e.g., Mathews & Brighenti 2003 and 
Bregman & Parriott 2009), as well as gas infall from the intergalac- 
tic medium. 

The hot gas in massive ETGs (^ lO^Afo) has been found 
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to radiate powerfully at X-ray wavelengths through thermal 
bremsstrahlung and yet it does not appear to be cooling as expected. 
The X-ray spectral energy distributions (SEDs) of massive ETGs 
demonstrate that hot gas {kT = 0.3-1 keV) typically dominates the 
0.5-2 keV emission (e.g., Boroson et al. 201 1). The temperatures 
and densities in the central regions imply relatively short radiative 
cooling times of » 10® yr (Mathews & Brighenti 2003). However, 
large quantities of cool (10‘*”^K) gas are not observed in the ETGs, 
which would be predicted by simple cooling flow models (see, e.g., 
Mathews & Brighenti 2003 for a review). By these observational 
arguments, it is necessary that some feedback mechanism (e.g., the 
AGN radio jets predicted by the models) keeps the gas hot and/or 
expels the cooled gas reservoirs. 

Direct observational evidence for the interaction of AGN ra- 
dio jets with hot gas has been obtained via X-ray and radio ob- 
servations of massive ETGs in the local universe (e.g., Boehringer 
et al. 1993; Birzan et al. 2004; Forman et al. 2005; Rafferty et al. 
2006). These observations have revealed relativistic radio outflows 
inflating large X-ray-emitting gas cavities with cool gas observed at 
the cavity rims. Measurements of X-ray cavity sizes and their sur- 
rounding gas densities and temperatures can give estimates of the 
mechanical energy input required by the radio jets to inflate the cav- 
ities against the pressure of the surrounding gas. The derived me- 
chanical energy and jet power are in the ranges lO^^-lO®^ ergs and 
10'*'’-10‘*® erg s~', respectively (McNamara et al. 2009, Nulsen et 
al 2007; astro-ph/0611136), sufficient to suppress gas cooling in 
the galaxy and impede star formation and cold gas SMBH accre- 
tion (Allen et al. 2006). This type of feedback has been referred to 
as ‘radio mode’ or ‘maintenance mode,’ since during this phase the 
accretion rate onto the central SMBH driving the AGN is low, and 
thus the majority of the gas heating is through radio AGN activity. 

These previous studies clearly indicate that heating by radio 
jets is an important process in galaxy evolution. Investigations of 
the influence of radio AGN on gas cooling in the general ETG pop- 
ulation (Best et al. 2005; hereafter BOS) found that more than 30% 
of the most massive (~ 5 x lO'^M©) galaxies host a radio AGN, 
which is consistent with radio AGN activity being eposidic with 
duty cycles of »10’-10® yr (Best et al. 2006). Evidence of such 
episodic radio luminous activity is also implied by the presence 
of multiple bright rims and shocks in the X-ray and radio images 
of individual ETGs (e.g., M87; Forman et al. 2005). Therefore the 
prevention of the cooling of large quantities of gas is thought to be 
maintained by a self-regulating AGN feedback loop. Cooling gas 
in the ETG centre initially provides a slow deposition of fuel for 
SMBH accretion, which in a radiatively-ineflicient accretion mode, 
leads to the production of a radio outburst. Surrounding cool gas 
is physically uplifted by the radio outbursts, which increases the 
gravitational potential energy of the gas or removes it from the sys- 
tem entirely (e.g. Giodini et al. 2010). As the gas further cools via 
X-ray emission, it falls back towards the SMBH where it can re- 
ignite a new cycle of accretion, thus completing the feedback loop 
(Best et al. 2006; McNamara & Nulsen 2007). The importance of 
the role of feedback from moderately radio luminous AGN is be- 
coming increasingly apparent, since the feedback energy can be 
directly diffused into the interstellar medium (Smolcic et al. 2009). 

A more complete understanding of the history of gas cool- 
ing and feedback heating in the massive ETG population requires 
direct X-ray and radio observations, respectively, of distant ETG 
populations covering a significant fraction of cosmic history. At 
present, such studies are difficult due to the very deep X-ray obser- 
vations required to detect the hot X-ray emitting gas in such distant 
populations (however, see, e.g., Ptak et al. 2007 and Tzanavaris 


& Georgantopoulos 2008 for some early work). Notably, Lehmer 
et al. (2007) utilised X-ray stacking techniques and the rj 250 ks 
Extended Chandra Deep Field South (E-CDF-S) and f»l Ms 
Chandra Deep Field-South CDF-S to constrain the evolution of 
hot gas cooling (via soft X-ray emission) in optically luminous 
{Lb ~ 10'® Z.g,a) ETGs over the redshift range of z « 0-0.7. This 
study showed that the mean X-ray power output from optically lu- 
minous ETGs at z ~ 0.7 is ~l-2 times that of similar ETGs in the 
local universe, suggesting the evolution of the hot gas cooling rate 
over the last k 6.3 Gyr is modest at best. Considering the relatively 
short gas cooling timescales for such ETGs ('--lO® yr), this study 
provided indirect evidence for the presence of a heating source. 
Lehmer etal. (2007) found rapid redshift evolution for X-ray lumi- 
nous AGNs in the optically luminous ETG population, which given 
the very modest evolution of the hot gas cooling, suggests that AGN 
feedback from the radiatively-efficient accreting SMBH population 
is unlikely to be the mechanism providing significant feedback to 
keep the gas hot over the last ;«6.3 Gyr. However, the mechanical 
feedback from radio AGNs, which is thought to be an important 
AGN feedback component, was not measured. 


In this paper, we improve upon the Lehmer et al. (2007) re- 
sults in the following key ways: (1) we utilise significantly im- 
proved spectroscopic and multiwavelength photometric data sets 
to select hot gas dominated optically luminous ETGs (via rest- 
frame colours, morphologies, and spectroscopic/photometric red- 
shifts) and sensitively identify AGN and star-formation activity in 
the population (see § 2 and 3); (2) we make use of a larger collec- 
tion of Chandra survey data (totaling a factor of «4 times the Chan- 
dra investment used by Lehmer et al. 2007) from the «2 Ms Chan- 
dra Deep Field-North (CDF-N; Alexander etal. 2003), the w4 Ms 
CDF-S (Xue et al. 201 1), and the «250 ks E-CDF-S (Lehmer etal. 
2005) (collectively the CDFs), which allows us to study the prop- 
erties of hot gas (e.g., luminosity and temperature) in optically lu- 
minous ETGs to z « 1 .2; and (3) we make use of new radio data 
from the VLA to measure the radio luminous AGN activity and the 
evolution of its duty cycle in the ETG population and provide di- 
rect constraints on the radio jet power available for feedback. The 
paper is organised as follows. In § 2, we define our initial work- 
ing sample and discuss the ancillary multiwavelength data used to 
identify non-passive ETG populations. In § 3, we use various se- 
lection criteria to identify passive ETGs and ETGs hosting radio 
AGNs. In § 4, we constrain the evolution of the X-ray emission 
from hot gas in our passive ETG sample using X-ray stacking tech- 
niques. In § 5, we discuss the level by which radio AGN can pro- 
vide heating to the hot gas in the ETG population. Finally, in § 6, 
we summarize our results. Throughout this paper, we make use of 
Galactic column densities of 7 Vh= 1 .6 x 10*® cm"^ for the CDF-N 
(Lockman 2004) and A'h= 8.8 x 10*® cm** for the E-CDF-S region 
(which also includes the CDF-S; Stark et al. 1992). In our X-ray 
aialyses, we make use of photometry from 5 bands: the full band 
(FB; 0.5-8 keV), soft band (SB; 0.5-2 keV), soft sub-band 1 (SBl; 
0.5-1 keV), soft sub-band 2 (SB2; 1-2 keV) and hard band (HB; 
2-8 keV). The following constants have been assumed, Qm = 0.3, 
Qa = 0.7 and //q = 70 km s**Mpc~' implying a lookback time 
of 8.4 Gyr at z = 1 .2. Throughout the paper, optical luminosity 
in the S-Band (Lb) is quoted in units of B-band solar luminosity 
(Lb,. ==5.2xl0®^ergs”'). 
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Figure 1. (a); Histogram of the entire sample of sources after imposing the selection criteria (i) to (iii) in § 2, showing the bimodal distribution of late-type 
and early-type galaxies. The dashed line shows the Bell et al. (2004) colour-cut, whereas the solid line shows the cut that we apply, slightly offset due to 
our different selection methods, using the redshift dependency taken from Bell et al. (2004) but varying the constant. The dotted and solid curves correspond 
to fitting a double gaussian to the sample in order to separate out the late and early-type populations, where the solid curve is the combination of the two 
gaussians. (b): Rest frame U-V colour versus absolute (Vega) F-band magnitude. The small filled grey circles represent the whole orginal sample of sources, 
again after imposing the selection criteria (i) to (iii) in § 2 (including both ETGs and LTGs). The open symbols represent the 393 sources in our final sample 
of ETGs and the dotted lines represent the colour cut applied using the minimum source redshift and maximum source redshift in our sample, with the line for 
the median redshift shown as a solid line (this cut was then applied to our model data later in Fig. 4). 


2 EARLY-TYPE GALAXY SAMPLE SELECTION 

The primary goals of this study are to constrain the potential heat- 
ing from AGN activity and the cooling of the hot gas in the opti- 
cally luminous (massive) ETG population over the redshift range 
z = 0.0-1. 2 (i.e., over the last 8.4 Gyr). To achieve these goals op- 
timally, we constructed samples of optically luminous ETGs in the 
most sensitive regions of the CDFs. 

We began our galaxy selection using master optical source 
catalogues in the CDF -N and E-CDF-S, which contain a collec- 
tion of IR-to-optical photometric data and good redshift estimates 
(either spectroscopic redshifts or photometric redshifts). The CDF- 
N master source catalogue consists of 48,858 optical sources de- 
tected across the entire CDF-N region (see Rafferty et al. 2011). 
This catalogue is based on optical sources detected in the Hawaii 
HDF-N optical and near-IR catalogue from Capak et al. (2004), 
and includes cross-matched photometry from GOODS-N through 
HST ACS and Spitzer IRAC catalogues (e.g., Giavalisco et al. 
2004), GALEX photometry,* and deep Aj-band imaging (Barger 
et al. 2008). In the E-CDF-S, we made use of a master catalogue of 
100,318 sources (see Rafferty et al. 2011). This catalogue is based 
on the MUSYC (Gawiser et al. 2006), COMBO-17 (Wolf et al. 
2004), and the GOODS-S (Grazian et al. 2006) optical surveys, and 
includes cross-matched photometry from MUSYC near-IR (Taylor 
et al. 2009), SIMPLE Spitzer IRAC (Damen etal. 2011), GALEX 
(see footnote 1), and GOODS-S deep [/-band photometry (Nonino 
et al. 2009). Our master catalogs are estimated to be complete 
to R< 26 (see section 2.1 of Xue et al. 2010). 

Whenever possible, we utilised secure spectroscopic redshifts, 
which were collected from a variety of sources in the literature and 


see http://galex.stsci.edu/GR4/. 


incorporated into the master source catalogues discussed above.^ 
When spectroscopic redshifts were not available, we made use 
of high-quality photometric redshifts, which were calculated by 
Rafferty et al. (2011) using an extensive library of spectral tem- 
plates (appropriate for galaxies, AGNs, hybrid galaxy and AGN 
sources, and stars), the optical-to-near-IR photometry discussed 
above, and the Zurich Extragalactic Bayesian Redshift Analyzer 
(ZEBRA; Feldmann et al. 2006). We compared these redshifts to 
the photometric redshift catalogue of Cardamone et al. (2010) 
finding a median difference of 0.01±0.16 between z = 0.0 — 0.8 
and O.OliO.32 between z= 0.8 - 1.2 in the two catalogues, thus 
providing additional evidence for the validity of these redshifts. 

Starting with the master catalogues of 149, 176 collective CDF 
sources, we imposed a series of selection criteria that led to the cre- 
ation of our optically luminous ETG catalogue that we use through- 
out this paper; the imposed selection criteria are summarized be- 
low: 

(i) We restricted ETG catalogue inclusion to sources with HST 
optical magnitudes of zg 5 o < 23 that were measured to be cosmo- 
logically distant (i.e., z > 0.05). The requirement of 2850 < 23 
ensures that the photometric redshifts of the remaining sources 
are of high quality ^ and provides a highly optically complete 
(see Fig. 2) sample of relatively bright optically luminous ETGs 
out to z ~ 1.2. Note that these photo-zs were computed using a 
redshift training procedure that implements spectroscopic redshifts. 
The true accuracy of the photometric redshifts is expected to be < 
6-7 times worse than those available for sources with spectroscopic 


^ For a comprehensive list of spectroscopic references, see Rafferty et al. 
( 2011 ). 

^ CDFN: median | zipec -Zphot I /(I ^Zspec) ~ 0.015, mean~0.032 and 
dispersionssO.090; E-CDF-S; median j z^pec -Zphot i /(• fZspec) « 0.007, 
meanwO.016 and dispersionwO.046, for zgso < 23 sources. 
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redshifts (see Luo et al. 2010 for details). This requirement further 
restricts our study to sources where HST imaging is available, thus 
allowing for further visual inspection of the optical morphologies 
to reasonably good precision (see criterion v below). This criterion 
restricted our working sample to 9732 CDF sources. 

(ii) We required that the sources are located within 6 ' of at least 
one of the six independent Chandra aimpoints in the CDFs (i.e., the 
«4 Ms CDF-S, four f«250 ks E-CDF-S pointings, and «2 Ms CDF- 
N). This criterion ensures that the galaxies are located in regions 
where the Chandra imaging is most sensitive and of highest quality 
(e.g., in these regions the Chandra point-spread function is small 
and relatively symmetric). Applying this additional restriction led 
to a working sample consisting of 6446 CDF sources. 

(iii) Using the redshift information available, we restricted our 
galaxy sample to include only sources with z = 0.05 — 1.2. The up- 
per redshift limit corresponds to the maximum distances to which 
we could obtain a complete sample of optically luminous ETGs 
that were relatively bright (2350 < 23) and contain useful morpho- 
logical information from HST imaging (see also, e.g., Haussler 
et al. 2007). Furthermore, this redshift upper limit for our survey 
allows us to detect the majority of X-ray luminous AGNs with 
Li-i kev ^ erg s”' located in the »2 Ms CDF-N and w4 Ms 
CDF-S surveys. This therefore defines the redshift baseline over 
which we can reliably measure hot gas emission through X-ray 
stacking without significant impact from undetected AGNs (see 
§ 4). As we will discuss below, when performing X-ray stacking 
analyses, we further exclude galaxies with z ^ 0.6 in the more 
shallow rs 250 ks E-CDF-S based on the same logic. For the mo- 
ment, however, our galaxy sample includes E-CDF-S sources with 
z « 0 . 6 - 1 . 2 , since we will later use these galaxies to constrain the 
radio AGN duty cycle in the ETG populations (see § 5.1). The 
imposed redshift limits led to the inclusion of 5734 galaxies with 
z= 0.05 -1.2. 

(iv) Since we are ultimately interested in measuring the hot gas 
X-ray emission from massive ETGs, we required that the galaxies 
that make up our sample have rest-frame 5-band luminosities in 
the range of Lb = 3-30 x 10 ’°Ib,g- As noted by O’Sullivan et al. 
(2001; see also Ellis & O’Sullivan 2006 and Boroson et al. 2011), 
such optically luminous ETGs in the local Universe have rela- 
tively massive dark matter halos, and are therefore observed to have 
0.5-2 keV emission dominated by hot interstellar gas {kT ^ 0.3- 
1 keV) with minimal contributions from other unrelated X-ray 
emitting sources (e.g., low-mass X-ray binaries; see Fig. 3b). This 
further restriction on including only optically luminous galaxies led 
to 243 1 galaxies. 

(v) To identify passive ETGs in our sample, we made use of the 
miiltiwavelength photometry and redshift information discussed 
above to measure rest-frame U — V colours, and we further used 
HST imaging to provide morphological information about our 
galaxies. As noted by Bell et al. (2004), the rest-frame U —V colour 
straddles the 4000 A break and provides a sensitive indication of 
mean stellar age. For our sample, we first required that all galaxies 
have rest-frame U — V colours redder than 

( 6 '-F)rea = 1.15-0.31z-0.08(iUr+22,4), (1) 

where My is the absolute F-band magnitude. Equation 1 (estab- 
lished to be valid out to z ^ 1 ) is based on Bell et al. (2004; see 
§ 5); however we have used a different constant term based on 
our analysis in Fig. la where we determine the red/blue galaxy bi- 
modal division by fitting a double gaussian to the distribution of 
U — V + 0.3 Iz + O.OSMk for our sample of galaxies after impos- 
ing the selection criteria (i) to (iii). In this exercise, we applied the 


redshift dependency from Bell et al. (2004) but shifted the con- 
stant (by 0.13) to fit to our sample, which is consistent with 
a typical colour scatter of < 0.2 mag for the red sequence colour- 
magnitude relation (see § 4 in Bell et al. 2004). We classified galax- 
ies lying below this divide as ‘blue cloud’ galaxies and those above 
as ‘red sequence’ galaxies. We (A.L.R.D. and B.D.L.) then visu- 
ally inspected all red-sequence galaxies using grayscale HST im- 
ages from the Z 850 band, and for the subset of sources located in 
the GOODS-N and GOODS-S footprints, we also inspected HST 
false-colour images based on 5435 , V^oe, and zgso observations. We 
strictly required the galaxies to have bulge dominant optical mor- 
phologies for ETG catalogue inclusion, and we rejected ETG can- 
didates that appeared to be possible edge-on spirals, which may 
simply be reddened by intrinsic galactic dust. Furthermore, we re- 
moved five sources which were very near the edges of the HST im- 
ages, where morphological classification was not possible. Apply- 
ing these morphological criteria led to our final sample of 393 opti- 
cally luminous ETGs. The basic properties of our parent sample 
are shown in Table 1. 

We note that out of the 393 optically luminous ETGs that make 
up our sample, 190 of these galaxies lie in the CDF-S or CDF-N at 
z = 0.05-1.2, or in the E-CDF-S at z = 0.05-0.6, which could po- 
tentially be used in X-ray stacking. The remaining 203 sources lie 
in the E-CDF-S at z = 0.6-1 .2. Of the 393 galaxies in our sample, 
163 have spectroscopic redshifts, and the remaining 230 sources 
have high-quality photometric redshifts from the Rafferty et al. 
(201 1) catalogue. Furthermore, 128 out of 190 sources potentially 
to be used in X-ray stacking have spectroscopic redshifts. 

Using both the photometric and spectroscopic redshifts we 
carry out a basic test of the environment of our sources by 
searching for neighbouring galaxies with an angular separa- 
tion of <500kpc from each of our 393 ETGs and within a red- 
shift difference of 0.09 and 0.046 in the CDF-N and E-CDF- 
S respectively (the dispersion in the photometric redshifts; see 
footnote 3). We first apply the cut in zgso optical magnitude of 
zg 5 o < 23 in order to ensure we are only using high quality pho- 
tometric redshifts and apply a further cut in absolute magni- 
tude of M 2850 < — 21. This results in a total galaxy sample of 
~ 3500 galaxies in the CDF-N and ~ 1000 in the E-CDF-S. We 
find a median of 10 ± 1 and 4 ±2 companions per ETG in the 
CDF-N and E-CDF-S respectively. We then check the number 
of neighbours we find for a random galaxy by searching within 
the comparison galaxy catalogues, and find a slightly lower me- 
dian number of companions within SOOkpc of 7 ± 2 and 3 ± 1 for 
the CDF-N and E-CDF-S respectively. Therefore, we find tenta- 
tive evidence that the massive ETGs are in richer than average 
environments (likely small groups). Since we are using the pho- 
tometric redshifts there is quite a large uncertainty, however, 
when using only the spectroscopic redshifts (giving us a much 
smaller and likely imcomplete sample) we do still find evidence 
for clustered environments. This is not unexpected since we are 
selecting massive ETGs. 


3 MULTIWAVELENGTH CHARACTEMSATIONS OF 
ETGS USING ANCILLARY DATA 

In this section, we make use of the extensive multiwavelength data 
available in the CDFs to identify both passive and non-passive (e.g., 
star-forming and AGN) ETGs. In the next section (§ 4), we will 
perform X-ray stacking analyses of the passive ETG population to 
measure directly the evolution of the mean hot gas emission. In the 
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Table 1. Master Catalogue, 


RA 

Dec 

z 

spec/phot? 

Z850 

Me/ 

Mb 

Mr 

Lb 

M. 

X-ray? 

1,4GHz? 

24/mi? 

(J2000) 

(J2000) 







log(Ls,0) 

log(Ms) 




(I) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

(8) 

(9) 

(10) 

(11) 

(12) 

(13) 

52.8483000 

-27.9371400 

0.816 

P 

21.52 

-21.79 

-22.06 

-22,83 

11.01 

11.16 

0 

0 

0 

52.8506205 

-27.9442900 

1.056 

P 

21,79 

-21.41 

-21.78 

-22.36 

10.91 

11.00 

0 

0 

0 

52.8527205 

-27.7069500 

0.526 

P 

20.62 

-20.49 

-20.83 

-21.65 

10,52 

10.84 

0 

0 

0 

52,8637605 

-27.6886300 

0,908 

P 

21.32 

-21.22 

-21.79 

-22.53 

10.91 

11,29 

0 

0 

0 

52.8672000 

-28.0023100 

0.727 

P 

21.78 

-20.64 

-20.96 

-21.55 

10.57 

10.73 

0 

0 

0 

52.8714105 

-28.0047900 

0,727 

P 

20.82 

-21.33 

-21.56 

-22.24 

10.81 

11,04 

0 

0 

0 

52.8717405 

-27.9800800 

0.771 

P 

21.81 

-21.25 

-21.55 

-22.15 

10,81 

10.88 

0 

0 

0 

52.8731805 

-28.0159400 

0.685 

P 

20.95 

-20.95 

-21.33 

-22.23 

10.72 

11.10 

0 

1 

0 

52.8741105 

-28.0181600 

0,727 

P 

21.64 

-20.81 

-20.99 

-21.67 

10.59 

10.80 

0 

0 

0 

52,8809595 

-27.7222000 

1,005 

P 

22.21 

-21.27 

-21.40 

-21.98 

10,75 

11,02 

0 

0 

0 


Notes: Columns (I)-(2): Optical J2000 coordinates. Column (3): Source redshift. Column (4): s=spectroscopic redshift, p=photometric redshift. Column 
(5): Zgso magnitude. Column (6): U-band magnitude. Column (7): B-band magnitude. Column (8): V-band magnitude. Column (9): Logarithmic S-band 
optical luminosity (log(ifl,;;)). Column (10): Logarithmic stellar mass derived from K-band magnimde (log(M©)). Column (1 1): Indicates whether the source 
was X-ray detected or not (0=not detected, 1 “detected). Column (12): Indicates whether the source was 1 ,4GHz radio detected or not (0=not detected, 
l=detected). Column (13): Indicates whether the source was 24/rm detected or not (0=not detected, l=detected). Table 1 is presented in its entirety (393 
sources) in the electronic version of the journal. Only a portion (first 10 sources) is shown here for guidance. 
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Figure 3. (a): The ratio of count-rates in the SB2 and SB 1 X-ray bands against redshift. The solid line at T^ff = 2.3 shows our adopted divide between hot-gas 
dominated sources and possible low mass X-ray binary (LMXB) and/or AGN dominated sources. The SED for a 1 .5 keV Raymond-Smith plasma (Raymond 
& Smith 1977) is shown as a dashed curve. Our stacked results (open triangles) are consistent with this SED. The errors on the band ratio were determined 
following the ‘numerical method’ described in § 1.7.3 of Lyons (1991). Band ratios are corrected for differential vignetting between the different bands using 
the appropriate exposure maps, (b): The rest-frame k-corrected SB (0.5-2 keV) X-ray luminosity (derived from the 0.5-1 keV flux to minimise contribution 
from LMXBs, and converted to 0.5-2keV luminosity using a 1.5 keV Raymond-Smith plasma SED) versus the faded, rest-frame S-band luminosity (where 
the evolution oflg is removed, using Faber etal. (2007) to parameterise the effect: Z,b,z=o = Sb x lo ® j[je jqIjj and the Icr dispersion (shaded) 

show the contribution to X-ray emission from hot gas (O’Sullivan et al. 2001; converted from bolometric to soft-band X-ray luminosity). The dashed line 
represents the expected contribution to X-ray emission from LMXBs (O'Sullivan et al. 2001); arrows represent the 3a upper limit X-ray luminosities. Those 
sources identified as hot-gas dominated from their band ratios are represented here as filled circles, and they lie well within 30 of the hot gas relation from 
O’Sullivan et al. (2001). The results from the stacking procedure in six redshift bins are shown as open triangles and these results are based on the SBl (0,5- 
1 keV) emission, convolved with the SED of a 1.5 keV Raymond-Smith plasma (Raymond & Smith 1977) to determine the corresponding SB (0.5-2 keV) 
emission. 


analyses below, we match our 393 optically luminous ETG optical 
source positions to those provided in multiwaveiength catalogues 
using closest-counterpart matching, which is a reasonable method 
provided that the optimum matching radius is carefully selected. 
We selected the optimum matching radius for each multiwave- 
length catalogue by first performing matching using a 30" match- 


ing radius and then observing the distribution of closest-counterpart 
matching offsets. For all catalogues discussed below (i.e. optical-x- 
ray, optical-radio, optical-infrared matching), we found the distri- 
bution of offsets to peak close to r»0", reach a minimum at sslf'5, 
and subsequently rise toward larger offsets due to spurious matches. 
A matching radius of lf'5 was therefore adopted as the optimum 
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Table 2. Hot gas dominated X-ray detected galaxies. 


RA 

Dec 

z 

SBl counts 

SB2 counts 

SB2/SB1 


LiaGHz 

log(L^.z=o) 

stack(y/n) 

(J2000) 

(J2000) 


(0,5-1 keV) 

(l-2keV) 


10“” ergs-* 

10”/4ly 

log(Lij.o) 


0) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

(8) 

(9) 

(10) 

03:32:09.706 

-27:42:48.110 

0.727 

74.367'|« 

125 

1 44+03! 

185.79l]5']f 

56,54±2,79 

10,90 

N 

03:32:28.734 

-27:46:20,298 

0.737 

61.891«« 


1 c'j-fo.io 

lO.nlj;’] 

59.63T2.81 

10.62 

Y 

03:32:34.342 

-27:43:50.092 

0.668 

40.2277fS 


i.z.z„o,29 

c iq+1.08 
:>.iy_0 92 

<5.39 

10.50 

Y 

03:32:38.786 

-27:44:48.923 

0.736 


<15.31 

<1.61 

1 .1+0.87 

13.10T1.42 

10.34 

Y 

03:32:41.406 

-27:47:17.185 

0.685 

7.971^1 

8.4li|;2 

1 A2+1-00 

1.06l»;« 

<5.72 

10,65 

Y 

03:32:44.088 

-27:45:41.461 

0.488 


1 1 1 

“•^^-4.62 

^ -^^-0.67 

0,49l»;i 

<2.56 

10.67 

Y 

03:32:46.536 

-27:57:13,104 

0.770 

21 1^9.69 

-8.50 

jQ5g 

1 19+0.67 
0.44 

68.68 :_23_o8 

<7,55 

10.95 

N 

03:32:46.949 

-27:39:02,916 

0,152 

24.771*;” 

<25.36 

<0.73 

1.32lgj^ 

0,65±0.04 

10.64 

Y 

03:32:52,066 

-27:44:25.044 

0.534 

12.67ip’ 

23 03 

1 1 -2+034 
‘•‘•^-0.54 

I7.88l|»] 

15.66±1.29 

10.87 

Y 

12:36:39.760 

62:15:47.832 

0,848 

14.42l«7 

<11.56 

<0,86 

Q 41+3.13 
*•^^-2.45 

<6.34 

10.49 

Y 

12:36:44.414 

62:11:33.347 

1.013 

12.111|| 

14.0015.” 

1.22+“]5 

1 2r+03§ 
*-’^0.32 

0 97+3-95 

876,37±28,74 

10.76 

Y 

12:36;52,895 

62; 14:44, 152 

0.321 

35.05+”] 

44.531*;” 

1 -24+0..^» 
-0.24 

6.41±0.31 

10,35 

Y 


Notes; Columns (1M2); Optical J2000 coordinates. Column (3): Source redshift. Column (4): (0.5-1 keV) net counts. Column (5): (1-2 keV) net counts. 
Column (6): (1-2 keV)/(0.5-l keV) count-rate ratio (SB2/SB1). Column (7): Rest-frame 0.5-2 keV luminosity (ergs s"') derived from SBl counts and 
Raymond-Smith plasma SED. Column (8): Radio luminosity (L\a ohz). Column (9): Logarithm of the S-band optical luminosity (logLs^e), Column (10): 
Indicates whether the source was included in our stacking analyses (Y/N). 



Redshift (z) 

Figure 2. The optical B-band luminosity of the ETG sample versus the 
redshift. The entire sample of 5734 sources after imposing the selection 
criteria (i) to (iii) in § 2 is represented by filled grey circles with the final 
ETG sample of 393 galaxies as open circles. The boxed areas show the bins 
from which sources were selected for X-ray stacking analysis, described in 
§4,1. The bins were selected to be evenly separated in co-moving volume. 


matching radius for all but the radio catalogues, for which the po- 
sitional errors are very small, therefore l"was more appropriate. 
Matches were visually inspected to further ensure they were sensi- 
ble, The number of spurious matches was determined for each data 
set analytically by calculating the ratio between the total area cov- 
ered by the parent sample sources, each with 1 "5 or l"matching ra- 
dius (Jt(1.5)* X 393 ~ 2778sq. arcsec or ti(l.O)^ x 393 1234sq. 

arcsec respectively) and the total area of the CDFs (within 6 ' of 
each pointing; 0.188 sq. degrees or 2436480 sq, arcsec). This ra- 
tio was then multiplied by the total number of sources in the mul- 
tiwavelength catalogues that lie within 6' of one of the Chandra 
aimpoints. 



Figure 4. Selecting passive galaxies to use in the X-ray stacking analysis 
based on their IRAC 5.8 fira/S /jm colours. The solid curve shows the me- 
dian 5,8 /m/8 /m colours of the 5000 model passive galaxy SEDs generated 
by hyperz with the 20 dispersion shown as dotted curves. In this figure, the 
most active galaxies (e.g. star-forming galaxies) have redder colours than 
those expected for passive ETGs. The error bars shown at z — 0.05 are 
representative error bars, plotted at the median 5.8/m/8/m colour for the 
following three cases: sources lying within ±20 of the passive line; out- 
liers with lower signal to noise with bluer colours (> 20 above the passive 
line); redder AGN and active star-forming galaxies (> 20 below the passive 
line). In our stacking analysis in § 4. 1 we test the effect of including/not in- 
cltiding the sources lying in regions which are > 2o away from the passive 
line but wliich have not otherwise been classified as active, however, there 
is very little difference in the final result aside from signal-to-noise. In our 
final stacking analysis presented in Fig. 7 these sources have therefore been 
included. 
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3.1 X-ray Properties of ETGs 

The ultradeep Chandra data in the CDFs provide a direct means 
for classifying X-ray detected ETGs as either hot gas dominated or 
likely AGNs. We used the published main catalogues for each of 
the CDFs, which consist of 503 sources in CDF-N (sk2 Ms; «0. 12 
deg^ survey, Alexander et al. 2003), 740 sources in the CDF-S 
(~4 Ms; »0.13 deg^ survey, Xue etal. 2011), and 762 sources 
in the E-CDF-S (four contiguous «250 ks Chandra observations 
that flank the CDF-S proper; jsO. 31 deg^, Lehmer et al. 2005). Us- 
ing our sample of 393 ETGs, the optical coordinates of the galax- 
ies were matched to the CDF X-ray catalogue positions using our 
adopted matching radius of 1 "5. When an ETG matched to a source 
in both the E-CDF-S and the CDF-S simultaneously (due to overlap 
between the E-CDF-S and CDF-S) we chose to use the data for the 
CDF-S, since these X-ray data are significantly deeper with smaller 
positional errors. In total, 55 X-ray matches were found once rep- 
etitions had been removed, including 1 1 in the CDF-N and 44 in 
the E-CDF-S region. The fraction of spurious matches in all the 
CDFs together was estimated to be =s2.8% (or ^2 expected spuri- 
ous matches). 

In Fig. 3a, we show the SB2/SBI count-rate ratio versus red- 
shift for the X-ray detected ETGs in our sample. The SB2/SB1 ra- 
tio provides an effective discriminator of the X-ray spectral shape 
in the SB, the energy regime where hot gas is expected to domi- 
nate. T 5 fpically, 0-2 AGNs have r^jf ~ 1 .8-2.3 (e.g., Alexander 
et al. 2005, Vignali et al. 2002, Reeves & Turner 2000). There- 
fore we took the upper limit of this range and conservatively clas- 
sified sources with SB2/SB1 ^ 1.7 (corresponding to Feg- 2.3) 
as sources having SB emission dominated by a hot gas compo- 
nent. Sources detected only in SB2 (i.e., having only a lower- 
limit on SB2/SB1), that had SB2/SB1 limits below our adopted 
cut were not classed as hot gas dominated sources. Sources with 
SB2/SB1 hardness ratio greater than this cut (i.e., SB2/SB1 > 1.7), 
have X-ray emission likely dominated by low mass X-ray binaries 
(LMXBs) or X-ray AGNs. However, by construction, our choice to 
study optically luminous ETGs will inherently minimise contribu- 
tions from LMXB-dominated systems and therefore AGNs are ex- 
pected to dominate the SB2/SB 1 > 1 .7 population (see below). Our 
SB2/SB1 criterion indicated 12 hot-gas dominated sources and 25 
likely AGNs (Fig. 3a). The SB2/SB1 ratios imply that a Raymond- 
Smith plasma (Raymond & Smith 1977) of kT—1.5keV is a good 
spectral model from which to convert count-rates to flux. In Table 2, 
we tabulate the properties of these X-ray detected ETGs. 

In Fig. 3b, we show the 0.5-2 keV luminosity (hereafter. Lx) 
versus Lb,q (see § 4 for details) for the ETGs in our sample. In 
order to minimise the contribution from LMXBs we calculated the 
rest-frame 0.5-2 keV luminosities Lx, SB based on the 0.5-1 keV 
SBl fluxes provided in the Chandra catalogues and convert them 
to 0.5-2 keV SB fluxes, applying a k-correction: 

Lx = 4rt4/xL (ergs^'), (2) 

where di is the luminosity distance in cm, /x is the 0.5-2 keV flux 
in units of erg cm“* s~*. The quantity k is the redshift-dependent 
k-correction. For sources that were characterised as hot gas domi- 
nated we used the observed 0.5-1 keV flux and a Raymond-Smith 
plasma SED (with kT\ ~ 1-5 keV; Raymond & Smith 1977; see 
Fig. 3a) to compute Lx- For sources that were identified as AGN 
dominant, we used a power-law SED (with F = 1.8) and the ob- 
serv'ed 0.5-2 keV flux to compute Lx- The solid line and shaded 
region shows the best-fit local relation and lo dispersion for hot 
gas dominated ETGs, and the dashed line shows the expected con- 


tribution from LMXBs (based on O’Sullivan etal. 2001 and typi- 
cally a factor of ~10 below the hot gas contribution). We note that 
nearly all ETGs without X-ray detections (plotted as upper limits) 
and the 12 ETGs with SB2/SBI band ratios consistent with be- 
ing hot gas dominated (highlighted with filled circles) also have 
LxJLb values similar to those observed for local hot gas domi- 
nated ETGs. The majority of the remaining X-ray detected sources 
with SB2/SB1 > 1.7 are expected to be AGNs. As Fig. 3b shows, 
these sources typically have large values of Lx! L b, again consis- 
tent with that expected from AGNs (see O’Sullivan et al. 2001; 
Ellis & O’Sullivan 2006). To further check for AGNs in our sam- 
ple we cross-matched our optical catalogue with spectroscopic data 
from Szokoly et al. (2004), Mignoli et al. (2005), Ravikumar et al. 
(2007), Boutsia et al. (2009) and Silverman et al. (2010) using a 
1 "5 radius in order to identify any sources with spectral features 
indicative of AGN, such as broad emission lines. We identified two 
potential broad line AGN in the E-CDF-S, both of which were X- 
ray detected and had already been flagged as likely AGN using our 
band ratio analysis (Fig. 3). 

3.2 5.8-24 /an Properties of ETGs 

In order to explore further whether the ETGs contained more subtle 
signatures of AGN or star-formation activity than provided by their 
X-ray and optical spectroscopic properties, we utilised Spitzer pho- 
tometry over the 5.8-24 pm range. We began by using Spitzer IRAC 
5.8 pmJi pm colours. Since AGN tend to be redder than galaxies 
in the mid-infrared, the 5.8 /an/8 pm colour can be used to iden- 
tify AGNs when the continuum is dominated by a rising power law 
component rather than a dropping stellar component (e.g., Stem 
et al. 2005). Similarly, the SEDs of powerful star-forming galax- 
ies, containing a large hot dust component, may exhibit this rise 
towards redder wavelengths and may also be identified by their 
5.8 /an/8 /an colour. In the E-CDF-S, we take the photometry in 
these channels from Damen et al. (201 1) and in the CDF-N, we take 
photometry from GOODS-N, cutting both catalogues at a signal-to- 
noise level of S/N< 5. 

We matched the positions from our sample of 393 ETGs to 
the Spitzer IRAC catalogues using a matching radius of 1/5 and 
found 384 matches. We estimated a spurious matching fraction of 
«3.1% (a;12 matches). In Fig. 4, we show the 5.8 /an/8 pm colour 
versus redshift for the 384 ETGs in our main sample. To determine 
the expected 5.8 pmj% pm colours for passive galaxies, we used 
the code hyperz (Bolzonella et al. 2000) with the SED library 
of Bruzual & Chariot (1993) to generate 5000 model galaxy SEDs 
based on a wide range of star formation histories and redshifts. For 
each galaxy, we adopted a formation redshift randomly selected to 
lie between the galaxy redshift (i.e., /age = 0) and 2 « 5. In Fig. 4, 
we plot the running median of the 5.8 /an/8 pm colours for the 
hyperz sample, after imposing the rest-frame U — V color crite- 
rion in equation 1 (the solid line in Fig. 4), and calculate the 2a dis- 
persion either side of the median (the dashed lines on Fig. 4). This 
curve was calculated by binning the data into bins of Ac = 0.05 
and computing the median and dispersion for each bin. Sources 
with very red IRAC colours lying below the passive line are likely 
AGN or star forming galaxies, and it can be seen that X-ray, ra- 
dio and 24 /an detected galaxies tend to lie below the solid line. 
In our stacking analyses (see § 4.1), we experimented with remov- 
ing the sources that lie outside of the 2o dispersion boundaries of 

^ http://data.spitzer.caltech.edu/popular/goods/ 
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Figure 5. Radio images of the three extended radio matches discovered by visual inspection of the radio images. The axes give the distance to the centre of 
the image in arcseconds. All of these extended sources have significantly higher 1.4GHz fluxes (1.2^.8mJy) than the other 16 radio-bright AGN in our ETG 
sample. 


our hyper-z normal galaxy envelope. However, we found no dif- 
ference in the general results since most of the sources exhibiting 
non-passive activity (either due to star formation or AGNs) had al- 
ready been identified by other indicators, and we therefore decided 
to include all of our galaxies in our subsequent stacking analyses 
(unless otherwise flagged as non-passive). As an additional test 
we experimented with the IRAC colour-colour diagnostic as in 
Stern et al. (2005), Fig. 1, however, we find that only three of our 
sources lie in their region of active sources, all three of which 
we have already flagged as active sources through our other di- 
agnostics. 

To identify additional ETGs in our sample that have signatures 
of star formation or AGN activity from dust emission, we cross- 
matched our ETG sample with Spitzer MIPS 24 pm catalogues. The 
E-CDF-S was observed with SpitzerlMlVSi as part of the FIDEL 
legacy program ^ (PI: Mark Dickinson; see § 2. 1.2-2. 1. 3 of Mag- 
nelli et al. 2009). We have used a catalogue of 20329 sources pro- 
duced by the DAOPHOT package in IRAF (see § 2.3 of Biggs etal. 
201 1). The MIPS 24 pm sensitivity over the E-CDF-S varies signif- 
icantly across the 30' x 30' field, with exposure times ranging from 
11,000-36,000 s. We make use of sources having signal-to-noise 
of at least 5a («30-70 /dy limits; see Magnelli etal. 2009). For 
the CDF-N, we made use of the publicly available GOODS Spitzer 
Legacy survey catalogues of 1 198 sources (PI: M. Dickinson). We 
utilised the 5a sample (flux limits of 70 /dy in the E-CDF-S and 
30 /dy in the CDF-N; Magnelli et al. 2009). Using a lf'5 match- 
ing radius we found a total of 20 matches to the 393 ETGs in our 
sample; three in the CDF-N and 17 in the E-CDF-S, with ~1.7 spu- 
rious matches expected. 24 pm provides a robust diagnostic of the 
presence of cold dust emission from the circumstellar envelopes of 
young embedded UV-luminous stars, characterised by a rising SED 
through the mid-infrared Muzerolle et al. (2004). Such systems are 
expected to contain significant X-ray contributions from popula- 
tions that are unrelated to hot gas, and we therefore classify these 
20 sources to be star-formation active systems. 


- http://irea.ipac.caltech.«iii/data/SPlTZER/FIDEL/ 


3.3 Radio Properties of ETGs 

To measure powerful radio emission produced by either radio-loud 
AGNs or star-formation activity, we cross-matched the optical co- 
ordinates of the parent sample with 1 .4 GHz VLA catalogues in 
the CDFs (using a 1" matching radius). We utilised the catalogue 
from Miller et al. (2008), but included additional sources at S/N>5 
(Miller et al. in preparation) for the E-CDF-S region, which con- 
tains 940 radio sources with S/N> 5 and reaches a 5o limiting flux 
density of 30 /Jy with a synthesised beam of 2,8" x 1.6". For the 
CDF-N, we utilised the catalogue from the Morrison et al. (2010) 
GOODS-N observations, which provides entries for 1227 discrete 
radio sources with S/N> 5 and 5a flux density limit of 20 /Jy at the 
field centre, with a 1.7" beam. In total, 24 radio detected counter- 
parts to the 393 ETGs were found (six in the CDF-N and 18 in the 
E-CDF-S) and 1 5 of these radio detected ETGs were also X-ray de- 
tected. The spurious matching fraction was estimated to be ajO.5% 
(fsO.l matches) and therefore negligible. Since the radio emission 
from radio luminous AGNs can be extended (e.g., FanarofF & Riley 
1974), the radio maps were carefully inspected by eye against the 
1" radius matching circles (overlaid at the locations of the parent 
sample positions) to verify the accuracy of the matches and iso- 
late extended sources. We identified three bright extended sources 
that were all identified using closest-counterpart matching; radio 
images of these sources have been provided in Fig. 5. We note 
that some of these individual sources have been well studied in 
the literature (e.g. J123644.4; Richards et al. 1998, J033238.8 and 
J033210.1; Kellermann et al. 2008). 

We calculated rest-frame 1 ,4 GHz monochromatic luminosi- 
ties for all radio detected sources following, 

I,,4OHz = 4rt4/,,4GHzl0'^®(l+z)““‘WHz-‘, (3) 

where /1.4 ghz the 1.4 GHz flux density (/dy) and a is the ra- 
dio spectral index for a power-law radio SED (i.e., Fy “c v““). 
We adopted a power-law spectral index of a = 0.85 (see Richards 
2000 for motivation). For normal galaxies without active radio 
AGNs, radio emission originates from HII regions and Type II 
and Ib supemovae, which produce synchrotron radiation from rela- 
tivistic electrons and free-free emission (Condon 1992). In passive 
ETGs, the contribution from these processes is unlikely to exceed 
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Table 3. Summary of source classifications. 


Classification 

No. of galaxies 

ETG sample 

393 

X-ray detected 

55 

Passive X-ray detected 

12 

Potential LMXB/X-ray AGN 

43 

Radio detected 

24 

Radio AGN 

19 

24/an detected 

20 

Star-forming galaxies 

15 

X-ray stacked galaxies (main) 

158 

X-ray stacked galaxies (faded) 

60 


(Ledlow 1997), The radio luminosity for all ra- 
dio detected sources in our sample was greater than 10^® W , 
which is expected given the flux limits of our survey. Therefore, de- 
tecting them at all suggests an excess of non-passive activity from 
either star- formation (SFR ^ 0,1 Mq yr“*) or AGN activity. 

To discriminate between star-formation and AGN activity in 
the radio-detected population, we use the well-known strong cor- 
relation between radio and far-infrared emission, which extends to 
cosmologically significant redshifts (at least z r; 1; Appleton et al, 
2004, and z> 2 using total infrared luminosity; Mao et al, 2011), 
For all our ETGs that are detected at both 24 /an and 1.4 GHz we 
measured the quantity /^24 = log(724/<m/ /i.4GHz) (Appleton et al. 
2004) (where aud /i. 4 ghz are observed fluxes). Radio-excess 
AGN can be identified by comparing their infrared emission to their 
radio emission, as their radio emission is significantly brighter than 
their infrared emission when compared to star-forming galaxies, 
which fit tightly along the far-infrared/radio correlation. Follow- 
ing Del Moro et al. (submitted) demonstrating the typical ^24 of 
radio-excess AGN based on starburst SEDs, we apply a selection 
of q 2 A < 0.5 to be indicative of radio AGN. This results in 19 of 
the 24 radio-detected galaxies being classified as radio AGN, with 
the remaining five radio-detected galaxies being classified as star- 
forming galaxies (as indicated in Table 4). In this exercise 10 ETGs 
with 24 /an detections but not radio detections were excluded from 
the final sample. We note that 16 of the 24 radio AGN were also 
X-ray detected. Table 4 shows the matched radio sources that are 
classified as radio AGN from the q 24 analysis, and which are used 
to estimate the AGN heating in § 5. We note that this approach 
identified all the sources with extended radio emission in Fig. 5 as 
radio AGN. 

In Table 3, we summarise the various source classifications 
described in § 3 for clarity. Of the original 393 galaxies in the ETC 
sample 190 of them can potentially be used in the X-ray stack- 
ing (see § 2). However, through various classification schemes we 
find that 32 of these are non-passive (potential X-ray AGN or star- 
forming galaxies) and are therefore excluded from the main stack- 
ing sample, thus leaving a sample of 158 passive galaxies which 
are suitable for X-ray stacking analyses. Of the 393 ETGs, 24 are 
radio detected and 20 of these are likely radio AGN while the other 
five have radio emission dominated by star formation. We classify 
a further 1 0 sources as likely star-forming galaxies, which have de- 
tections only in 24 /an and not radio, and lower limits of q 24 > 0.5. 
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Figure 6. Ratio of observed flux in 24 /an and radio (1.4GHz) versus red- 
shift, Sources which lie below the solid line (Del Moro et al, in prep) are 
taken to be radio AGN, 


4 COSMIC HISTORY OF X-RAY EMISSION FROM 
MASSIVE ETGS 

Approximately 90% of the ETGs in our passive sample are unde- 
tected in the X-ray (338 galaxies). Therefore to measure the hot gas 
emission from the whole population it is necessary to implement 
X-ray stacking techniques. This investigation focuses on measur- 
ing the cooling of the hot gas in ETGs, which dominates emission 
at soft X-ray energies (0,5-2 keV), as opposed to LMXBs, which 
dominate emission in the hard band (2-8 keV). Stacking analyses 
were therefore carried out in the soft bands (i.e., SBl and SB2), 
which we expect to be dominated by hot-gas emission and to have 
minimal contributions from LMXBs. In Fig. 2, we plot six redshift 
intervals of galaxies with Lg = (3 — 30) x 10'®Z,g^3 where we per- 
formed stacking analyses for each subsample (solid boundaries). 
The redshift divisions were chosen to encompass roughly equal in- 
tervals of comoving volume, and the larger redshift interval spacing 
beyond z « 0.6 is the result of excluding from our stacking analy- 
ses sources that were within the «250 ks exposure of the E-CDF-S 
(see §2 for details). 

We note that previous studies (e.g.. Bell etal. 2004; Faber 
etal. 2007) have shown that, from z = 1 - 0, the B-band lumi- 
nosity of typical massive ETGs fades by ~1 mag. To estimate 
the mean X-ray luminosity evolution for an ETG population with 
similar z = 0 B-band luminosities (Lb.z=o), we thus constructed 
six “faded” redshift-divided subsamples of ETGs with Ba, 2 =o = 
(3 -30) X 10’®ij,3. We calculated Lb.z=o following the Faber 
etal. (2007) prescription: Lb,z=o==Lb x 10"® '*x 1 - 23 xz these 
faded luminosities only 60 of the 158 stacking sources lay within 
the allowed range of optical luminosities. In total, we stacked 12 
subsamples of ETGs (six main and six faded) with both the main 
and faded samples having the same divides in redshift but with 158 
sources in the total main sample and only 60 in the total faded sam- 
ple. 
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Table 4. Radio-bright ETGs. 


RA 

(J2000) 

(1) 

Dec 

(J2000) 

(2) 

z 

(3) 

/ l .4 GHz 

(/dy) 

(4) 

log(L|.4C/ft) 

(lO^^WHz”') 

(5) 

(/dy) 

(6) 

924 

(7) 

Lysa 

(10'*" erg s”'cm”^) 
(8) 

log(La) 

(log(La,0)) 

(9) 

Extended 

(Y/N) 

(10) 

Note 

(AGN/SF) 

(H) 

03:31:29.563 

-28:00:57.384 

0.685 

46.3±8.8 

0.9+0.2 

<70 

<0.18 


10.7 

N 

A 

03:31:32.210 

-27:43:08.076 

0,956 

68.7±7.3 

2.9±0.3 

85.5±6.2 

0.09 

<5.1 

10.9 

N 

A 

03:31:39.041 

-27:53:00.096 

0.220 

61,2±7.0 

0.08±0,01 

85.5±5.3 

0.15 


10.6 

N 

A 

03:31:40,044 

-27:36:47.628 

0.685 

208.1L15.4 

4.0±0.3 

<70 

<-0.47 


11.2 

N 

A 

03:31:45.895 

-27:45:38.772 

0.727 

42.7L6.8 

0.9±0.2 

<70 

<0.22 


11.3 

N 

A 

03:31:57,782 

-27:42:08.676 

0,665 

97,2±6.5 

1.7±0.1 

73.7±3,1 

-0.12 

8.4 

10.8 

N 

A 

03:32:09.706 

-27:42:48.110 

0.727 

257,2±12.7 

5.6±0.3 

<70 

<-0.57 

18.6 

11,3 

N 

A 

03:32:10,137 

-27:59:38.220 

0,685 

1165.0±36.0 

22.2±0.7 

<70 

<-1.22 


11.2 

Y 

A 

03:32:19.305 

-27:52:19,330 

1.096 

39,1 ±6.2 

2.3±0.4 

<70 

<0.25 


11.3 

N 

A 

03:32:28.734 

-27:46:20.298 

0.737 

263,3±12.4 

6.0±0.3 

<70 

<-0.58 

1.0 

11.0 

N 

A 

03:32:28.817 

-27:43:55,646 

0.214 

4814.0±103.0 

6.3±0.1 

<70 

<-1.84 

0.04 

10,6 

Y 

A 

03:32:38,786 

-27:44:48.923 

0.736 

58.0±6.3 

1.3±0.1 

247.8±2.4 

0,63 

0.2 

10.7 

N 

S 

03:32:39.485 

-27:53:01.648 

0.686 

107.0±6,2 

2.04±0,I2 

<70 

<-0,18 

<0,2 

11,0 

N 

A 

03:32:46.949 

-27:39:02.916 

0,152 

105.5L7.0 

0.060±0.004 

106.2±2.7 

0.003 

0.1 

10.7 

N 

A 

03:32:48.177 

-27:52:56.608 

0.668 

32.8L6.2 

0.6±0.1 

152.4±2.3 

0.67 

4.9 

11,2 

N 

S 

03:32:52.066 

-27:44:25.044 

0.534 

148,3±12.2 

1.6±0.1 

<70 

<-0.326 

1,79 

11,1 

N 

A 

03:33:05,671 

-27:52:14.268 

0.521 

55,7±6.8 

0.6±0 ,l 

181,5±7.6 

0.51 

<2,4 

11.3 

N 

S 

03:33:15.427 

-27;45:24,012 

0,727 

63.5±6.9 

1.4±0.2 

368.3±2.7 

0.76 


10.6 

N 

S 

12:36:01.813 

62:11:26.659 

0,913 

99.2L5.5 

3.8±0.2 

<30 

<-0.52 


10.9 

N 

A 

12:36:08,137 

62:10:36.136 

0.679 

213.1±7.9 

4,0±0.2 

2300.0±12.9 

1.03 

0.2 

10,7 

N 

S 

12:36:17.098 

62:10:11.554 

0.846 

65.3±8.3 

2.1±0.3 

88.2±6.4 

0.13 

7.95 

10.7 

N 

A 

12:36:22,705 

62:09:46,313 

0,748 

45,7±5,1 

1.1±0.1 

<30 

<-0.18 


10,6 

N 

A 

12:36:44,414 

62:11:33.347 

1.013 

1805.1 ±59.2 

87.6±2.9 

<30 

<-1.78 

0.9 

11.3 

Y 

A 

12:36:52.895 

62:14:44,152 

0.321 

198.3±9.6 

0.64±0.03 

<30 

<-0.82 

0,1 

10.5 

N 

A 


NOTES; Columns (l)-(2): Optical J2000 coordinates. Column (3): Redshift, Column (4): Radio (l,4GHz) flux density (/dy). Column (5): Radio (1,4 GHz) 
luminosity (loglO^^W Hz”'). Column (6): 24 /m flux density (ijjy). Column (7): qit ratio, (log/24^„//l.4C//r). Column (8): 0.5-2 keV flux 
(10“*^ ergs s”' cm”^) derived from SBl (0.5-1 keV) counts and converted using the 1.5 keV Raymond-Smith plasma. Column (9): B-band luminosity 
(login, :;)• Column (10): Indicates whether there is extended emission (Y=yes, N=no). Column (1 1): Note on classification: S = radio detected ETGs for 
which the radio emission is likely dominated by star formation (as implied by their value), A = radio emission dominated by an AGN. 


4.1 X-ray Stacking Technique 

Our stacking procedure, summarized below, makes use of images, 
background maps, and exposure maps that were constructed by 
Alexander et al. (2003) for the ^2 Ms CDF-N, Lehmer et al. (2005) 
for the «250 ks E-CDF-S, and Xue etal. (2011) for the =e4 Ms 
CDF-S. 

We chose to use circular apertures of constant radii to extract 
on-source counts. We chose to extract X-ray counts (source plus 
background) from a 1 "5 radius circular aperture centered on the 
locations of sources that were within 6' of any of the six Chan- 
dra aimpoints (the ^2 Ms CDF-N, the r;4 Ms CDF-S, and the 
four «250 ks pointings in the E-CDF-S). These choices of source 
inclusion radius and extraction aperture radius were previously 
found to optimise the stacked signal (see, e.g., Lehmer et al. 2005b, 
2007) and are therefore implemented here. For each source, we 
used our source extraction aperture to extract source plus back- 
ground counts Si from images and exposure times f; from expo- 
sure maps. For each stacked sample, total source plus background 
counts were computed as X = X,-A' and exposure times were com- 
puted as r = 0.03 h (the 0,03 factor comes from the fact that r,- 
is the sum of exposure map values over ss30 pixels). 

Background and exposure maps were then used to measure the 
background counts and exposures for each source. For this exercise, 
we used a 15" radius circular aperture centred on the location of 
each source to extract local background counts h;,iocai exposure 


times t/, local- The on-source background counts 6,-, were estimated 
following bj = local h/h',iocai- Total stacked background counts 
were then obtained through the summation B = bi- 

For each stacked sample, any galaxy that was classified as a 
normal ETG (see §3) was stacked. The stacking procedure was 
carried out with three different samples: (a) a sample including 
all radio AGN, passive X-ray detected sources (10 galaxies) and 
passive X-ray undetected sources; (b) a sample including only X- 
ray undetected galaxies and (c) a sample including passive X-ray 
undetected galaxies and radio AGN but excluding passive X-ray 
detected galaxies. However, we found that the inclusion of both X- 
ray detected normal galaxies and radio AGN in the stacking did not 
significantly change our results, implying that most X-ray luminous 
AGN had been successfully excluded from the sample via direct X- 
ray detection and classification. Therefore all radio AGN and X-ray 
detected normal galaxies (thus all passive galaxies) were included 
in all of the stacks resulting in a final sample of 158 and 60 passive 
ETGs to be stacked in the main and faded samples, respectively. 

For each stacked sample, we measured the signal-to-noise ra- 
tio (S/N = (X - B)/VB). For a significant detection, we required 
that S/M > 3, and for such stacks, we measured net counts as 
N = S — B. When a stacked sample was not detected, we placed 

3a upper limits on the net counts i.e. N < 3 + w|. where a^ 

and aB are the bootstrapped errors on the total and background 
counts respectively. The error on the net counts was determined 
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by applying a bootstrapping method. For each stacked sample con- 
taining n galaxies, we randomly drew n sources from the sample 
(allowing for multiple draws of the same source) and restacked the 
scrambled sample to measure net counts. This exercise was per- 
formed 1000 times for each stacked sample, thus giving a sense of 
the variance of the population. The count-rates O for each stacked 
sample were determined as <I> = t,N/T, where ^ is a mean aperture 
correction. Since many of the sources with relatively large PSFs (at 
~3-6' off axis) had aperture radii that did not encompass the whole 
PSF, it was necessary to factor in a correction (^,) for each of the 
stacked sources. The average correction factor used was computed 
as ^ s X(^/h7Z/h' where and t,- are the correction factors and 
exposure times measured for each individual source. The stacked 
count-rates O were then converted to fluxes using the SED for a 1.5 
keV Raymond-Smith plasma (Raymond & Smith 1977; see Fig. 3a 
for motivation and see Table 5 for kTx values used in each stack 
sample). Errors on the count-rate to flux conversion were calculated 
by propagating errors on the mean SB2/SB 1 ratio, which can be 
used as a proxy for temperature in our Raymond-Smith SED. The 
errors on the luminosity were determined by propagating the boot- 
strapped errors on the source counts and the errors on the countrate- 
to-flux conversion factor based on the errors in the SED tempera- 
ture described above. When calculating luminosities, the luminos- 
ity distance is calculated using the mean redshift in each bin (since 
as Fig. 2 demonstrates, the redshifts are quite evenly distributed in 
each stacking bin). 

4.2 X-ray Stacking Results 

The results of stacking the X-ray data of the sample, are sum- 
marised in Table 5. In Figures 3a and 3b, we have overlaid our 
stacking results, demonstrating that our galaxies have SB2/SB1 
count-rate ratios consistent with our adopted 1.5 keV Raymond- 
Smith plasma SED (which we use to convert count-rates to fluxes 
and luminosities) and Lx versus Lb,z=o values consistent with local 
hot gas dominated ETGs. In Figures 7a and 7b, we display the rest- 
frame 0.5-2 keV luminosity (computed following equation 2) per 
5-band luminosity LxILb versus redshift for our main and faded 
samples, respectively {circles). The 10 individually X-ray detected, 
hot gas dominated ETGs are shown as triangles in Figures 7a and 
7b (only eight are shown on Fig. 7b as only eight fulfilled the con- 
straints of the stacking for the faded sample). Only 10 out of the 
12 galaxies shown in Table 2 are included in the stacking as, for 
stacking, we add the limitation that all galaxies in the E-CDF-S 
must have z < 0.6. As would be expected from an X-ray selected 
subset, these sources generally have higher values of Lx/Lb- For 
comparison, we have also plotted the mean Lx/Lb values obtained 
by Lehmer etal. (2007) for ETGs with Lb > 10'° Lb.q (squares). 

To constrain evolution to z = 0, we take the Boroson 
et al. (2011) sample of 30 nearby ETGs and select only 
those 14 with Lb = (3-30) x10'®Ib, 3- We convert their 0.3- 
8 keV luminosities to 0.5-2 keV luminosities using our adopted 
1.5 keV Raymond-Smith plasma SED and find a mean value of 
logLx/LB = 29.7 ± 0.2 (crosses in Figs. 7a and 7b). The com- 
bination of the Boroson etal. (2011) mean Lx/Lb and our stack- 
ing results indicates that there is little apparent evolution in Lx/Lb 
for these optically luminous ETGs. However, a Spearman’s p test 
reveals that the quantity Lx/Lb is correlated with z at the 92% 
and 96% probability level for the main and faded samples re- 
spectively. To constrain the allowable redshift evolution of Lx/Lb, 
we fit a simple two parameter model to the data logix/^B = 
A +Slog(l +z) and find best-fit values of [A,B] = [29.73 ±0.14, 


1.19±0.68] and [29.71 ±0.14, 1.07 ±0.68] for the main and faded 
samples, respectively. These values indicate mild evolution in the 
X-ray activity of luminous ETGs and are consistent with those of 
Lehmer et al. (2007). Using this model, we find that at z=1.2, ETGs 
are ~ 2.4 ± 0.9 times or ~ 2.2 ± 0.9 times (for the main and faded 
sample respectively) more X-ray luminous (per unit Eg) than at 
z=0, which suggests only modest evolution. Our best-fit relations 
have been highlighted in Figures 7a and 7b as dashed curves. 

Since the stacked X-ray properties (i.e., SB2/SB1 band ratio 
and Lx versus Lb) are consistent with those expected from hot gas 
dominated ETGs, with little expected contributions from LMXBs, 
we can use the X-ray luminosity versus redshift diagram for our 
stacked samples as a direct tracer of the hot gas cooling history 
for massive ETGs with = (3 - 30) x 10'°Ls,©. The observed 
mild decline in X-ray luminosity per unit 5-band luminosity and 
roughly constant X-ray gas temperature for massive ETGs over 
the last »;8.4 Gyr of cosmic history suggest that, on average, the 
gas is being kept hot. We expect that many complex processes are 
contributing to the evolution of the gas including radiative cool- 
ing, periodic AGN heating and outflows, replenishment from stellar 
winds and supemovae, interactions and sloshing, and intergalactic 
medium and poor group inflow (e.g. Tabor & Binney 1993, Best 
et al. 2006, Faber & Gallagher 1976 and Brighenti & Mathews 
1999). The detailed influences that each of these processes has 
on the gas are difficult to quantify, particularly without a strong 
idea of the environment in which each galaxy resides. However, we 
know that most of our sources reside within small groups and clus- 
ters, therefore, processes such as intergalactic medium and poor 
group infall may be important. One of the goals of this paper is to 
test whether AGN feedback from mechanical feedback can provide 
enough energy to keep the gas hot and counter the observed cooling 
over the long baseline of cosmic time spanned by our observations. 

In the next section, we discuss the viability of AGN feedback 
heating of the gas by directly measuring the history of radio AGN 
events in our galaxy population and computing the mechanical en- 
ergy available from these events. 

5 DISCUSSION 

5.1 The Hot Gas Cooling and Mechanical Heating Energy 

Budgets 

The above X-ray stacking results indicate that the X-ray power out- 
put from hot gas in the massive ETG population remains well reg- 
ulated across a large fraction of cosmic history (since 2 ~ 1), To 
determine whether the heating from AGNs is sufficient to keep the 
gas hot, we estimated the mechanical power input from AGNs and 
the radiative cooling power from the hot gas. As discussed in § 4.2, 
the history of gas cooling power can be directly inferred from our 
X-ray stacking results; the gas cooling power, can be expressed as 

Ecool =CboiIx ~CboiiBl0-'^(l (4) 

where A = 29.73 ± 0.14 and 5 = 1.19 ± 0.68 were computed in 
§4.2, Lb » 6.3 X 10*® Lb.q is the mean value of Lb, and Cboi ~ 
1.8 is the bolometric correction for a hot gas SED with 1.5 keV 
temperature (see §4.2 above). In Fig. 8, we plot the mean cooling 
history (filled circles and dashed curve for stacked values and best- 
fit model, respectively), since z « 1 .2. 

To estimate the energy input from radio AGNs over the last 
fsS.4 Gyr of cosmic history, we began by measuring the radio AGN 
fraction as a function of radio luminosity (a proxy for mechanical 
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Table 5. X-ray stacking properties, 


Exposure time 




Net counts (S-B) 

SBl 

SB2 

S/N 

fjsi 

kTx 

La-.xs 

^B,mean 


^mean ^iot 

SBl 

SB2 

(Ms) 

(Ms) 

SBl SB2 

log(erg s“’cm"^) 

(keV) 

log(ergs“') 

\og{LB,(o) 

log(ergs”'L^‘.) 

(1) (2) 

(3) 

(4) 

(5) 

(6) 

(7) (8) 

(9) 

(10) 

(11) 

(12) 

(13) 


General Sample 


0.278 

19 

80,7±42,2 

125.3±66.3 

12.2 

12.3 

29.3 

36.4 

-16.1 


40,7±0,2 

10.7 

29.9±0.2 

0.519 

19 

38,8±20,3 

55.0±28.1 

25,7 

26.3 

9.9 

11.1 

-16.8 


40.6±0,2 

10.8 

29.8±0.2 

0.630 

33 

58.2±40.5 

97.3±50,2 

85.2 

85.5 

7,7 

10.2 

-17.0 

40.5±0.3 

10.7 

29.8±0.2 

0.790 

43 

112.6±59.0 

178.8±93.5 

111,6 

111.6 

12.7 

16,5 

-16.9 

1 5+0.3 

40,9±0.2 

10.8 

30.2±0.2 

0.973 

32 

35.4±12.9 

43.5±17,8 

71.7 

71,6 

5.0 

5.1 

-17,1 

1 

40,9±0.2 

10.8 

30.1±0.1 

1.124 

12 

6.1±4,8 

9.4±6.8 

34.7 

34,7 

1.3 

1.5 

<-17.2 

<40.9 

10.9 

<30.1 

Faded Sample 

0.284 

13 

29.3±23.2 

32,0±24.3 

6,2 

6.4 

15,8 

13.3 

-16.5 

1 

^•^-0.1 

40.3±0.3 

10,7 

29,7±0.3 

0.518 

11 

35,1±17,6 

51.9±24.5 

17.6 

18.2 

10.9 

13.0 

-16.8 

1 'I'l+O-OS 

40,5±0,2 

10.7 

29.9tli 

0,632 

13 

49.4±34.9 

79.4±48.4 

37.3 

37.7 

9.8 

12.6 

-17.1 

1 5+^-2 

40.4±0.3 

10.6 

29.8±0.2 

0.767 

15 

105.8±54.9 

151.0±85,2 

41,7 

41.7 

19.4 

22.3 

-16.9 

1 

^••^-0.4 

40.8±0,2 

10.6 

30,2±0.3 

0.980 

5 

16.5±8.5 

18.4±9,6 

13.4 

13.3 

5.4 

4.8 

-17.4 

1 2+0-5 

l.+ _0,4 

40.6±0.2 

10.7 

30.0±0.3 

1.119 

3 

2.6±3.1 

2.8±1.4 

9.6 

9,6 

1,0 

0.8 

<-17,5 

<40.7 

10.6 

<30.0 


Notes: Column (1): Mean redshift of bin. Column (2): Number of sources in stacking bin. Columns (3)-(4): Net counts for SBl and SB2. Columns (5)-(6): 
Exposure times in Ms for SBl and SB2. Columns (7)~{8): Stacked signal-to-noise ratio for SBl and SB2. Column (9): Logarithm of the stacked SBl flux 
(ergs s"* cm~^). Column (10); X-ray temperature (keV). Column (1 1): Logarithm of the stacked 0.5-2 keV luminosity Lx (ergs s""') measured using the 
SBl flux and an assumed Raymond-Smith plasma SED with 1.5 keV temperature. Column (12): Logarithm of the S-band luminosity (for the faded sample, 
we list Lb,z-o)- Column (13) Logarithm of the 0.5-2 keV to B-band luminosity ratio (for the faded sample, we list ix/f-B, 2 =o). 


M 30.5) 
Sh 


O 


X-ray detected sources 
Stacked result 
Boroson et al. 2011, hot gas 
Lehmer et al. 2007 


(a) 


0.0 


Faded X-ray detected sources 
Faded stacked result A 

Boroson et al. 2011, hot gas 
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Figure 7. (a): The evolution of the soft X-ray properties of hot-gas dominated ETGs in terms of B-band optical luminosity, including passive X-ray detected 
galaxies, passive X-ray undetected galaxies and radio AGN. All galaxies with individual detections in the SB band are plotted as triangles. The circles represent 
the stacked results in six redshift bins with only an upper limit in the highest redshift bin. There is very little evolution ofLysBi^^B in these hot gas dominated 
galaxies, in particular given the very short expected cooling time of hot gas. To minimise the contribution from LMXBs we stack in SBl (0.5-lkeV) but we 
use the SED of a 1 ,5keV plasma to convert this to SB (0.5-2keV) in order to more easily compare the data to previous studies, z = 0 points are shown for the 
hot gas contribution as determined from a subset of the ETG sample of Boroson et al. (2011). The simple two-parameter fit for logLx/Ln = A i-Blog(l +s) 
is shown as a dotted curve on both panels (a) and (b), with best-fit values of [T.B] = [29.73 ±0.14, 1,19 ±0.68] and [29.71 ±0.14, 1,07 ±0.68] for the main 
and faded samples, respectively, (b): Represents the faded sample, showing that the fading of Lb with cosmic time is a fundamental effect and acts to slightly 
decrease the observed Lx,sb/Lb ratio with increasing redshift. The errors on the luminosity are derived from propagating the bootstrapped errors on the net 
counts and systematic errors on the conversion between count-rate and flux using the l.SkeV Raymond-Smith plasma SED. 
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heating) and redshift. By making the assumption that all galaxies 
will go through multiple AGN active phases, we can use the radio- 
luminosity and redshift dependent AGN fraction as a proxy for the 
typical AGN duty cycle history for galaxies in our sample. 

To establish a baseline local (z « 0) measurement of the ETG 
radio AGN fraction, we used the BOS sample of radio-loud AGN 
from the SDSS survey, which included both early-type and late- 
type galaxies. For the sake of comparing these data with our ETG 
sample, we selected galaxies in the BOS sample with elliptical-like 
concentration indices C > 2.6 (Strateva et al. 2001). The concen- 
tration index is defined as C = rgo/rso, where rgo and r^o are radii 
containing 90% and 50% of the optical light respectively. By ap- 
plying the flux density limit of 5 mJy, we limit the BOS sample to 
a lower radio luminosity limit of «10^^ W Hz“*, which corre- 
sponds to a maximum redshift of 2 = 0. 1. To measure the AGN 
fraction for galaxies in the distant Universe, we used the sam- 
ple of distant ETGs presented in this paper. Using Z,i .4ghz = 
10^^ W Hz“*, the luminosity limit used for the BOS data, we 
determined that the corresponding CDF-N and E-CDF-S ra- 
dio flux limits (see §3.3) allow us to study similar AGNs out to 
z « 1 and 0.85, respectively. We calculated the AGN fraction 
for both the local BOS local galaxies and our distant galaxies in 
three bins of radio luminosity (in even logarithmic luminosity 
intervals) in the range of L\a gHz ~ (1 — 100) x 10^^ W Hz^', 
each bin with a different allowed redshift range due to the flux 
limits. These bins in luminosity and redshift then result In a 
total of 2642 elliptical galaxies containing a radio AGN from 
the BOS sample and 13 radio AGN from our sample. The AGN 
fraction was computed in each bin, for both the BOS sample 
and our sample, by taking the total number of radio AGN in 
a particular luminosity range and dividing it by the total num- 
ber of galaxies in which an AGN with a luminosity lying within 
that range could have been detected if present. We estimate la 
errors on the AGN fractions following Gehrels (1986). By com- 
paring the radio AGN fraction at the mean redshift z = Zmcan 
of our distant galaxy sample with that of the BOS z w 0 sample, 
we can estimate the evolution of the duty cycle of the AGN out- 
bursts in each radio luminosity bin. The time-dependent radio 
AGN fraction for each bin of radio-luminosity was computed 
following 



Figure 8. Ratio of mechanical power to 1.4 GHz radio luminosity 
(l-mech/t-i 4 GHz) versus Z.i.4 qhz for the 21 giant elliptical (gE) galaxies 
studied by Cavagnolo etal. (2010; filled circles). The Cavagnolo etal, 
(2010) relation, which includes gE galaxies and radio galaxies at the 
centers of cooiing clusters (from Birzan ct al 2008), has been shown as 
a dashed line. Our best-fit relation for gE only galaxies, presented in 
equation 6, is shown as a solid line. 


are expected to be gE and group central galaxies, we made use 
of Lmech and 1^1.4 GHz values for the sample of 21 gEs from Cav- 
agnolo et al. (2010) to derive the Z.mech~^i 4 GHz correlation for 
these sources. Figure 8 shows the 21 gEs from Cavagnolo et al. 
(2010). We find that the best-fit relation from Cavagnolo et al. 
(2010) {dashed line in Fig. 8), which includes radio galaxies at 
the centers of galaxies, overpredicts the ImechAi. 4 GHz ratios 
for AGNs with Ii. 4 ghz ~ 10^^ W Hz~’. Using these data, we 
derived the following relation, which is applicable to gE galax- 
ies: 

( 6 ) 


^ r ^ /agn,z=z™.u-/agn,z«0,, _ ^ ,,, Our best-fit relation is plotted in Figure 8 as a solid line. 

/4GA'IM.i,4 ghz,i) — t +/ agn,zk 0,)) (6) Using equations (5) and (6), we then estimate the average 

mechanical feedback power per galaxy over the last r; 8.4 Gyr 
where A/, is the difference in the mean lookback time between pf cogjjjic history considering all radio AGNs in the range of 
z == ^mean.i and z = 0 (i.e., our AGN fraction and that of BOS) in ^ ^ (1-100) x 10^^ W Hz"' via the following summation; 

a particular bin of mean radio luminosity 4 ghz./- 

Several studies have now shown that the radio power out- ^heating = '^fAGN{l<L\A GHz,/) ^mech.i (7) 

put Lj 4 GHz from AGNs within nearby giant elliptical (gE) and ' 

cluster central galaxies correlates with the inferred mechanical In Fig. 9, we show the mean heating luminosity and I 0 errors 

power imech th®* i* needed to inflate the cavities within hot X- {solid curve with shaded envelope) derived following equation 7. 

ray halos (e.g., Birzan et al. 2004; Birzan et al. 2008; Cavagnolo From Fig. 9, tve see that on average there appears to be more 

et al. 2010; O’Sullivan et al. 2011). Until recently, these rela- than sufficient input mechanical energy from radio AGN events 

tions have been calibrated using the cores of cooling clusters, to balance the hot gas radiative cooling. From the five stacked 

and may not be appropriate for the massive early-type galaxies bins where we obtain X-ray detections, we estimate on average 

studied here. Cavagnolo et al. (2010) have added a sample of 21 iheatingAcooi ~ 2.4 Iq This result is broadly in agreement with 

gE galaxies and have shown that, as long as the radio structures that found for local elliptical galaxies of comparable mass where 

are confined to the hot X-ray emitting gas region, gE galaxies the mechanical power has been measured using X-ray cavities (e.g., 

provide a natural extension to the imech~f"i.4 GHz correlation at Nulsen et al. 2007). Nulsen et al. (2007) estimate that the total cav- 

low £mech- However, as noted by Cavagnolo et al. (2010), gE ity heating can be anywhere between 0.25 and 3 times the total gas 

galaxies and FRI sources in group environments (e.g., Croston cooling if 1 pV of heating is assumed per cavity; however, the en- 

et al. 2008) tend to have Lnicch/ii.4 GHz ratios much lower than thalpy of the cavity and therefore the total heating may be much 

the correlation derived including clusters. Since our galaxies higher. Stott et al. (submitted) find a trend in groups and clusters 


©2011 RAS, MNRAS 000, ??-16 




14 Danielson et al. 


indicating the ratio of intra-cluster medium (ICM) AGN heating 
(from the brightest cluster galaxy) to ICM cooling increases with 
decreasing halo mass. For halo masses < 5 x IO'^Mq the heating 
can exceed the cooling. We have found that most of our galaxies 
are likely to live in small group environments. Therefore extrapo- 
lating this relation to the expected halo masses of the galaxies in our 
sample («[1-10] x 10‘2 m© for Lb » [3-30] x 10'° ig,©; Vale & 
Ostriker 2004) would similarly imply that the mechanical heating 
would likely exceed the radiative cooling. 

We note that our heating calculation is based on duty cycle 
histories derived primarily from the <20 distant radio AGNs in our 
sample and is based on the assumption that each galaxy will have 
many radio outbursts that span the full range of radio luminosities 
studied here. We therefore expect these calculations will have sig- 
nificant uncertainties that we cannot determine. In the next section, 
we estimate the global ETG hot gas cooling and radio AGN heating 
power as a function of redshift. 


5.2 Cosmic Evolution of Global Heating and Cooling Density 

Using a large sample of radio-loud AGNs in the local Universe, 
Best etal. (2006) computed the radio-luminosity and black-hole 
mass dependent AGN fraction of nearby galaxies. Their data 
show that the population-averaged mechanical power (probed by 
1.4 GHz power) produced by these AGN events increases with 
black-hole mass (and also S-band luminosity) and balances well the 
radiative power output ifom X-ray cooling of the hot gas (see their 
Fig. 2). Their analyses further revealed that relatively low luminos- 
ity radio AGN (logii ,4 ghz/(W Hz“*) « 22-25) are likely to pro- 
vide the majority of the mechanical feedback power for the popula- 
tion as a whole. They estimated that in the local universe, the mean 
mechanical power output density from mechanical heating from ra- 
dio AGNs with Ii ,4 GHz > 1Q22w Hz^‘ is ^4 x 10^' W Mpc-^. 

Due to the relatively small number of radio AGNs found in our 
survey, it is not feasible to calculate the evolution of the radio and 
mechanical luminosity density of the Universe. However, the evo- 
lution of the radio AGN luminosity function has recently been mea- 
sured out to z Ki 1.3 using the VLA-COSMOS survey (Schinnerer 
et al. 2007) to relatively faint luminosity levels (L\ a qhz 10^’- 
1 0^^ W Hz^ * ; Smolcic et al. 2009). By converting radio luminosity 
into mechanical luminosity, Smolcic et al. (2009) integrated their 
luminosity functions to determine the estimated mechanical power 
density of the Universe out to z « 1.3. In Fig. 10, we show the 
expected mechanical feedback power density evolution, based 
on the Smoliic et al. (2009) radio luminosity function and equa- 
tion 6 for both pure-luminosity density evolution (L» =« [1-fzp) 
and pure-density evolution (<1>, « [i 1), the best-fit parame- 
terisations for the evolution of the 1.4 GHz luminosity function. 

As shown in § 4, our X-ray stacking measurements can be de- 
scribed on average as ^cool.mod ^ Qjol 10^(1 Tz)®ig (with Cbol — 
1.8, .4 = 29.73 ±0. 14, B = 1. 19 ±0.68, and Lb « 6.3 x 10‘° Lb,s) 
for the ETG population. Using this scaling relation and the ob- 
served evolution of the ETG 5-band luminosity function (pg from 
Faber etal. (2007; see their Table 4), we can compute the expected 
volume-averaged cooling luminosity density. In this exercise, we 
assumed a constant intrinsic scatter of a 1 dex for the Lx/Lb 
ratio (Boroson etal. 2011) and transformed the 5-band luminos- 
ity function into a X-ray gas cooling luminosity function using the 
following transformation: 

(px(logix,z) = f q>j?(logLg,z)5(IogIx|logiB)rflogig 



0.0 0.2 0.4 0.6 0.8 1.0 1.2 

Redshift(2) 


Figure 9. {top panel) Mean radiative cooling power Icooi and mean me- 
chanical heating power iheating versus redshift. The filled circles and !o 
error bars show the bolometrically corrected estimates of Lcooi (see Table 5) 
and the solid curve shows our best-fit model. The tong-dashed curve with 
shaded region represents our best estimate of the mean heating luminosity 
as presented in equation 7. These measurements show that for the early-type 
galaxies in our sample, there is more than enough energy available from ra- 
dio AGN heating to keep the gas from cooling, {bottom panel) Ratio of 
mean heating to radiative cooling luminosity versus redshift. We find that 
the average heating power is rsl.4-2.6 times larger than the radiative 
hot gas cooling power. 


5(logix|log5B) 




(logTcool.mod ^OgLx)^ 
’2o2 


( 8 ) 


The total redshift-dependent cooling density fljooijz) of the Uni- 
verse can therefore be computed following 

4:ooi(z)=/ <l?x{'to%Lx,z)Lxd\ogLx. (9) 

Jo 

In Fig. 10, we show the resulting GcooiC^) versus redshift {solid 
curve). Our analyses show that the estimated mechanical power 
provided by radio AGN activity is a factor of «1.5-3.5 times 
larger than the radiative cooling power (see bottom panel of 
Fig. 10), and the shape of the heating and cooling histories appear 
to be in good overall agreement. For comparison, we plot the mean 
values of Lheating/icooi as measured in § 5.1 and Figure 9 {filled 
circles), which are in agreement with the global heating-to-cooling 
estimates obtained here. 

The combination of the approaches for estimating global and 
mean galaxy heating and cooling taken here and in § 5.1, respec- 
tively, indicate that mechanical heating exceeds that of the radiative 
gas cooling for early-type galaxies with Lb « (3-30) x 10*® Ta,©. 
These computations are based on the assumption that the radio 
luminosity provides a direct proxy for mechanical power, which 
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Redshift (z) 

Figure 10. (top pane!) Cosmic history of mechanical heating Cheating and 
radiative cooling ncooi of hot gas since z w 1.1. The mechanical luminos- 
ity density versus redshift for radio AGNs, as computed by SraolSid et al. 
(2009), for pure-luminosity evolution and pure density evolution are indi- 
cated as dashed and dotted curves, respectively. The shaded envelope pro- 
vides the estimated uncertainties between the i 1.4 qhz and correla- 
tion. The cooling history of hot gas in ETGs is indicated as a solid curve 
(see §5.2 for details), (bottom panel) Ratio of mechanical heating to radia- 
tive cooling versus redshift. For comparison, we have plotted our estimates 
of the mean heating-to-cooling luminosity ratios (iheating/icool as provided 
in the bottom panel of Fig. 9 (filled circles). This relation shows that there 
is plenty of mechanical energy potential available to keep the gas hot over 
the last s:8,4 Gyr. 


scales for our early-type galaxies in the same way as that measured 
for local gEs (i.e., based on data from Cavagnolo et al. 2010). 
Indeed, some studies have suggested that AGN heating in less mas- 
sive systems, like those studied here, may have different heating 
cycles and mechanical efficiencies (Gaspari et al. 2011). Future 
studies that characterize how the radio and mechanical power are 
related in galaxies like those studied here would be needed to ex- 
clude the possibility that the excess of mechanical power compared 
with cooling power (as observed here) is due to the calibration. 


6 SUMMARY AND FUTURE WORK 

The X-ray and multiwavelength properties of a sample of 393 mas- 
sive ETGs in the Chandra Deep Field surv'eys have been studied, 
in order to constrain the radiative cooling and mechanical feedback 
heating history of hot gas in these galaxies. We detected 55 of the 
galaxies in our sample in the X-ray bandpass, and using the X-ray 
and multi wavelength properties of these sources, we find that 12 of 
these systems are likely to be dominated by X-ray emission from 
hot gas. To measure the evolution of the average ETG X-fay power 


output, and thus hot gas cooling, we stacked the 0.5-1 keV emis- 
sion of the X-ray undetected and detected “normal” galaxy popula- 
tion in redshift bins. 

We find that the average rest-frame 0.5-2 keV luminosity 
per unit S-band luminosity (Lx/Lb) has changed very little since 
z « 1.2 and is consistent with (1 evolution. This sug- 

gests that the population average hot gas power output is well reg- 
ulated over timescales of ss8 Gyr; much longer than the typical 
cooling timescale of the hot gas (»0. 1-1 Gyr). We hypothesize that 
mechanical heating from radio luminous AGNs in these galaxies 
is likely to play a significant role in keeping the gas hot, and we 
compare the implied gas cooling from our stacking analyses with 
radio-AGN based estimates of the heating. 

We find that if local relations between radio luminosity and 
mechanical power hold at high redshifts, then the observed radio- 
luminosity dependent AGN duty cycle suggests that there would 
be more than sufficent (factor 1. 4-2.6 times) mechanical energy 
needed to counter the inferred cooling energy loss. Similarly, we 
find that the evolution of the mechanical power density of the 
Universe from radio AGNs increases only mildly with redshift and 
remains a factor of «1. 5-3.5 times higher than the radiative hot 
gas cooling power density of ETGs in the Universe, These results 
are concordant with previous lower redshift studies (e.g Best 
et al. 2005 and Lehmer et al. 2007) and with theoretical feedback 
models such as Churazov et al. (2005), Croton et al. (2006) Bower 
et al. (2008) and Bower et al. (2006) where feedback from radio 
AGN maintains the balance between heating and cooling rates of 
hot interstellar gas in massive ETGs. 

Understanding the evolution of both the X-ray and radio prop- 
erties of optically luminous ETGs could be improved in future work 
by: (1) gaining a better understanding of the environments of these 
sources in order to better constrain the likely contribution of e.g. 
gas infall to the evolution of the X-ray properties and probe the 
influence of environment on the balance between gas cooling and 
heating in galaxies. This could be achieved by measuring spectro- 
scopic redshifts for the whole sample, as there is currently signif- 
icant uncertainty in photometric redshifts. (2) Conducing deeper 
X-ray observations could provide stronger constraints on the evo- 
lution of the hot gas. We focus our observations in the soft band 
in which the background is lowest, therefore doubling the Chandra 
exposure time to 8Ms could provide a factor of 1 .4-1 ,6 improve- 
ment in the sensitivity resulting in the faintest detectable sources 
having soft band fluxes of 6.0 x 10"'*erg cm^^ s"’ (Xue et al. 
2011) and improved statistics on the average X-ray emission from 
the ETG population. 
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